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ABSTRACT 

We present millimeter molecular- line complemented by optical observations, along with a reanalysis 
of archival centimeter H I and continuum data, to infer the global dynamics and determine where 
dense molecular gas and massive stars preferentially form in the circumnuclear starburst ring of the 
barred-spiral galaxy NGC 7552. We find diffuse molecular gas in a pair of dust lanes each running 
along the large-scale galactic bar, as well as in the circumnuclear starburst ring. We do not detect 
dense molecular gas in the dust lanes, but find such gas concentrated in two knots where the dust 
lanes make contact with the circumnuclear starburst ring. When convolved to the same angular 
resolution as the images in dense gas, the radio continuum emission of the circumnuclear starburst 
ring also exhibits two knots, each lying downstream of an adjacent knot in dense gas. The results agree 
qualitatively with the idea that massive stars form from dense gas at the contact points, where diffuse 
gas is channeled into the ring along the dust lanes, and later explode as supernovae downstream of 
the contact points. Based on the inferred rotation curve, however, the propagation time between the 
respective pairs of dense gas and centimeter continuum knots is about an order of magnitude shorter 
than the lifetimes of OB stars. We discuss possible reasons of this discrepancy, and conclude that 
either the initial mass function is top-heavy or massive stars in the ring do not form exclusively at 
the contact points where dense molecular gas is concentrated. 

Subject headings: galaxies: individual (NGC 7552) — galaxies: ISM — galaxies: kinematics and 
dynamics — galaxies: starburst 



1. INTRODUCTION 



Over half (e.g., >7Q%: lEskridge et al.l[2QQQl ) of ah spi- 
ral galaxies in the Local Universe exhibit large-scale 
bars. Such barred galaxies exhibit circumnuclear star- 
bursts more frequently than non-barr ed spiral galaxies 
(jJogee et al.l [20051 : lEllison et al.l[20TTh . The circumnu- 
clear starburst usually takes the form of a ring ce ntered 
on the nucleus (e.g., review by 'Buta Combesl [^96|), 
and has a star- formation rate (SFR) that can make up a 
significant fr action, if not dominate, the total SFR of the 
galaxy ( Joge e et al.ll20Q5l ). Is gas transported from the 
surrounding regions towards the center of the galaxy to 
fuel the circumnuclear starburst, and if so how? If bars 
transport gas inwards, then why we do not see circum- 
nuclear starbursts in virtually all barred spiral galaxies? 
Stated in a different way, why is star formation is strongly 
promoted in the circumnuclear region of some galaxies? 
In this work, we address the last question through ob- 
servations of the barred spiral galaxy NGC 7552, which 
exhibits a circumnuclear starburst ring with a SFR of 
10-15 Mq yr~^ that comprises at least half of the total 
SFR of this galaxy. 

Theoretical studies suggest that, in galactic disks. 
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gravitational torques induced by bar potentials drive gas 
inflows. In this scenario, molecular gas experiences two 
principal shocks: (i) at the two leading edges of the 
large-scale bar, where dust lanes (indicating, presum- 
ably, the pileup of molecular gas) are commonly seen 
(e.g., Sheth et al. 2002); (ii) at the two locations where 
the inner terminal portions of the bar make contact with 
the circumnuclear starburst ring, where twin peaks in 
mole cular gas as traced by CO are commonly observed 
fe.g..'Ken nev et al.lll992l iSakamoto et al.l ll999^ During 
the passage of molecular gas in the galactic disk through 
the bar, angular momentum is transferred from the gas 
to the bar through the aforementioned shocks, resulting 
in a net inflow of gas along the dust lanes. The inflowing 
gas is eventually halted at an inner Lindblad resonance 
(ILR; or, if there are two ILRs, then between the inner 
ILR and outer ILR), which has the form of a ring, where 
the bar terminates and within which the net torque ex- 
erted by the bar is zero (|Buta fc Combesiri996l ). Because 
the formation of a circumnuclear gas ring requires a bar 
that is strong enough (i.e., to induce a sufficiently strong 
non-axisymmetric potential) to promote adequate gas in- 
fall, the mere presence of a bar is not sufficient for the 
presence of a circumnu clear gas ring. Simulations (e.g., 
iRegan fc TeubenI 12004), as well as observati ons (e.g., 
ISakamoto et al.l Il999l : iLaurikainen et al.1 l2007f ) , suggest 
that galaxies with long and massive (morphological clas- 
sification SBa) bars are more effective at channeling gas 
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towards their centers than those with weak bars (SBc), 
and hence preferentiahy host circumnuclear gas rings. 

Although molecular gas has been observed to accumu- 
late along the leading edges of the bar and at the ILR 
in agreement with model predictions, there is a strik- 
ing difference between the observed properties of the 
molecular gas and the star-formation activity in these 
two regions. The molecular gas in (at the leading edge 
of) the bar has been found to be much more diffuse 
on average (e.g. Downes et al. I [19961 : [Hiittemeister et al.' 
^000; Schinnerer et al. 2002), and the star- formation ac- 
tivity in the bar also much less pronounced, than that in 
the circumnuclear starburst ring. Based on observations 
of NGC 6951, a strongly barred spiral galaxy that ex- 
hibits a circumnuclear starburst ring and which in many 
other respects ar e remarkably simil ar to NGC 7552 as 
explained below, Kohno et al. ("1999") described how dif- 
fuse molecular gas may collect in the circumnuclear star- 
burst ring of this galaxy and the circumstances under 
which this gas subsequently becomes sufficiently dense 
to form stars. They found that, just like in a number 
of other barred galaxies having circumnuclear starburst 
rings, diffuse molecular gas as traced in ^^CO is concen- 
trated in two knots on opposite sides of the ring where 
a pair of dust lanes, each running along the (leading 
edge of the) bar on opposite sides of the nucleus, make 
contact with the ring, demonstrating that the twin CO 
peaks likely represent molecular gas channeled inwards 
along the dust lanes. For the first time, they showed 
that dense molecular gas as traced in HCN is not cospa- 
tial with the CO knots, but is instead concentrated in 
two k nots located ju st downstream of the respective CO 
knots. iKohno et al.l (jl999) argued that diffuse molecular 
gas, channeled inwards along xi orbits (closely traced by 
the dust lanes), is too turbulent to become dense where 
they intersect with X2 orbits (traced by the circumnu- 
clear starburst ring), as the gas at these locations (de- 
lineated by the CO knots) experience relatively strong 
shocks. Instead, they proposed that the molecular gas 
only becomes dense further downstream where the afore- 
mentioned turbulence has partially dissipated, leading to 
the collapse of ^as clouds vi a gravitational (Toomre-type) 
instabilities ()Toomrelll964f ) and subsequent star forma- 
tion. 

Because intercomparisons can provide insights not pos- 
sible from studies of a single object alone, here we test 
the idea proposed by iKohno et al.l ([1999) on another 
galaxy. NGC 7552 is a member of the Southern Grus 
Quartet at a distance of 22.3 Mpc (jKennicutt et al.ll2003L 
so that V = 107 pc). Observations of the entire group 
in atomic hydrogen (H I) gas have revealed an exten- 
sion on the south-western side of NGC 7552 indicative 
of tidal interactions with one or more group members 
(^Dahlem 2005), although no such distortion is visible 
in the optical disk suggesting that only the outskirt of 
NGC 7552 is strongly perturbed. The gross optical mor- 
phology of NGC 7552 — circumnuclear starburst ring, 
large-scale bar with dust lanes along the leading edges, 
and two prominent spiral arms emerging from the ends 
of the bar as can be seen in Figure [T] — is remarkably 
similar to that of NGC 6951, including even the orien- 
tation of the abovementioned features in the sky. The 
large-scale stellar disk of NGC 7552 is close to face-on, 
with an axial ratio between the optical major and minor 



axes indicating an inclination of about 25° (for compar- 
ison, also based on axial ratios from optical images, the 
LEDA extragala ctic database quotes a n average inclina- 
tion of 23°, and iFeinstein et all () 19901 ) infer an inclina- 
tion of 28°). The nuclear activity of NGC 7552 is clas- 
sified as a Low-Ionization Nu clear Emission-line Regio n 
(LINER) with H II reg ions (jPurret Bergeronlll987[ ). 
iClaussen fc Sahail ()1992l ) find that the molecular gas in 
NGC 7752 as traced in ^^CO (1-0) is strongly concen- 
trated in t he n u cleus a nd the bar. 

Forbes et al.l ()1994a[ ) have imaged the circumnuclear 
starburst ring of NGC 7552 in the radio continuum at 
6 cm and 3 cm, providing the least-obscured view of the 
ring so far available. As explained in greater detail later, 
the ring is nearly circular, and exhibits five prominent 
knots. The radio spectrum of each knot falls between 
6 cm and 3 cm, indicating nonthermal synchrotron emis- 
sion prod uced, presum ably, as a result of supernova ex- 
plosion s. iForbes et al.l (1994b) also imaged the ring in 
Ha and lSchinnerer et"ai.i (JJ97) in Br7, where both these 
lines trace gas ionized by recently- formed massive stars. 

Our work presented here involves multi-line radio in- 
terferometric complemented by optical observations of 
NGC 7552, along with a reanalysis of archival data for 
this galaxy in H I and centimeter continuum. In ^we 
provide details on the observations performed and the 
data reanalyzed, and in ^ the results obtained. In g] 
we derive the physical properties of the molecular gas, 
and in ^the global kinematics of NGC 7552. Then, in 
^ we examine the implications of our results and anal- 
yses for the formation of dense molecular gas and stars 
in the circumnuclear starburst ring of NGC 7552, and 
pr opose an imp ortant revision to the picture advanced 
by Kohno et alj (j 19991) . based on observations of the cir- 
cumnuclear starburst ring in NGC 6591, that is better 
able to explain the observed properties of the circumnu- 
clear starburst ring in NGC 7552 as well as more recent 
observations of dense gas in the circumnuclear starburst 
ring of NGC 6951. Finally, in ^ we provide a concise 
summary of our work. 

2. OBSERVATIONS AND DATA REDUCTION 

Our millimeter molecular-line observations toward the 
center of NGC 7552 focussed on four bright tracers of 
molecular hydrogen gas, ^^CO, ^^CO, HCN, and HCO+. 
^^CO traces relatively diffuse (~10^~^ cm~^) molecular 
gas characteristic of the bulk of giant molecular clouds 
(GMCs) in our galaxy, whereas HCN and HCO+ traces 
relatively dense (~10^~^ cm~^) molecular gas character- 
istic of dense condensations in GMCs where stars form; 
as we will show in §4.2[ in our observations ^^CO traces 
molecular gas at intermediate densities. The spatial 
distribution of ^^CO therefore indicates locations where 
atomic hydrogen gas has presumably been compressed to 
form molecular hydrogen gas, whereas the spatial distri- 
bution of HCN and HCO+ indicates locations where the 
molecular gas has attained densities sufficiently high to 
form stars. These molecular-line observations also pro- 
vide information on the gas kinematics: as described in 
§5.H we used the ^^CO image to infer a rotation curve 
for the central region of NGC 7552. The H I image of 
NGC 7552 provides information on the gas kinematics 
beyond the central region imaged in molecular gas, and 
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as described in ^5.11 was used to infer a rotation curve 
for the galaxy between the innermost to the outermost 
detectable regions of the H I disk. The rotation curve 
thus inferred was used to derive the predicted locations 
of dynamical resonances in the galaxy, as well as to com- 
pute dynamical timescales in the circumnuclear starburst 
ring. We compared the predicted locations of dynamical 
resonances to individual features in NGC 7552 from our 
wide- and narrow-band optical observations, thus pro- 
viding an assessment of the reliability of the inferred 
rotation curve. The narrow-band Ha image from our 
optical observations also provided estimates of the SFRs 
in different parts of the galaxy. Finally, from the radio 
continuum image of the circumnuclear starburst ring, we 
inferred locations where large concentrations of massive 
stars recently met their demise as super novae. 

2.1. Optical observations: Gemini South 

The optical observations of NGC 7552 reported here 
are the product of a winning entry in the 2011 Aus- 
tralian Gemini School Astronomy Contest, for which one 
of us (S. Ryder) was a co-investigator on the contest pro- 
posal. Images of the galaxy in g\ r\ and i' filters were 
obtained with the Gemin i Mult i- Object Spectrograph 
fGMOS: iHQok et al.l[2QQl attached to the Gemini South 
Telescope as part of the program GS-2011A-Q-12 (PI: 
C. Onken). The filters are similar to the filters used by 
the Sloan Digital Sky Survey (SDSS). We also obtained 
an Ha image with a narrow-band Ha filter centered on 
656 nm with a 7 nm wide bandpass. The observations 
were done in seeing between 0% — 0'!8 under photometric 
conditions on 2011 July 23 UT. Four exposures per fil- 
ter were obtained with integration time between 100 and 
300 s, with each exposure offset by 10 arcsec spatially to 
allow filling-in of the inter-CCD gaps in CMOS. On-chip 
binning yielded a pixel scale of 0'.'146 pix~^. 

The data were reduced and combined using the gem- 
ini package (VI. 10) within IRAF0. A master bias frame 
(constructed by averaging with 3-cr clipping a series of 
bias frames) was subtracted from all raw images. Im- 
ages of the twilight sky were used to flatfield the images, 
then the dithered galaxy images were registered and av- 
eraged together with the imcoadd task to eliminate the 
inter- CCD gaps, bad pixels, and cosmic rays. 

2.2. Radio observations: Australia Telescope Compact 

Array 

2.2.1. HI 

We re-created the 21-cm neutral atomic hydrogen (H I) 
maps of Dahlem (2005), who kindly provided us with cal- 
ibrated data, using robust weighting to provide a higher 
angular resolution than the maps he published. The pa- 
rameters of the H I observation are summarized in Table 
[H The H I maps presented here have a synthesized beam 
of 20'.'0 X 20'.'0 at full- width half-maximu m (FWH M) (by 
comparison, the H I maps pubhshed by iDahlemI ((2005) 
have a synthesized beam of 34Vl x 30V0). In our chan- 
nel maps, the root-mean-square (rms) noise level in each 
13 km s~^ channel is 1.1 mJy beam~^. The integrated 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



H I intensity and intensity-weighted mean H I velocity 
maps, as well as all other such maps presented in this 
work, were constructed using the task MOMNT in the 
AIPS package; this task allows us to select or discard 
pixels in each channel map based on whether their corre- 
sponding intensities when smoothed over a selected num- 
ber of channels exceed a specified threshold. The H I 
disk of NGC 7552 is much smaller than the size of the 
primary beam, and hence no primary beam correction 
was applied to the data. 

2.2.2. Radio Continuum 

We re-processed the radio continuum data taken by 
Forbes et al. (1994a) at 3 cm and 6 cm, and re-created 
their maps to permit our own analyses as well as pro- 
vide digital versions for use in making overlays. As 
summa rized in Table [TJ the observations of iForbes et al.l 
([1994^ ) were conducted in one 12-hour synthesis on 11- 
12 February 1993 using the Australia Telescope Compact 
Array (ATCA) in its 6C configuration. This configu- 
ration provides baseline lengths ranging from 153 m to 
6000 m, making it possible to image structures (with uni- 
form brightnesses) as large as about 285^.^0. The primary 
fiux calibrator used was PKS B 1934-638, for which we 
adopted a fiux density of 5.83 Jy at 6 cm (4786 MHz) and 
2.84 Jy at 3 cm (8640 MHz). The quasar PKS B2311-452 
was used for complex gain calibration. 

We calibrated the data in the standard manner using 
the MIRIAD package. To improve the dynamic range of 
our maps, we self-calibrated the datasets at both 6 cm 
and 3 cm from the maps obtained at the respective wave- 
lengths. The final map at 6 cm has a synthesized beam 
of 2^0 X 1V3 at a position angle (PA) of 0°, and an rms 
noise fiuctuation of 72 /iJy beam~^. The final map at 
3 cm has a synthesized beam of x V!l and an rms 
noise level of 66 /iJy beam~^. Both the angular resolu- 
tions and noise levels of our map s are, of cour s e, sim i- 
lar to those in the maps shown bv lForbes et al.l (|1994a[ ). 
Given the small angular extent of the circumnuclear radio 
continuum emission from NGC 7552, no primary beam 
correction was applied to the data. 

2.2.3. HCNandHCO'^ 

We observed NGC 7552 in HCN (1-0) and HCO+ (1-0) 
with the ATCA on three separate occasions between May 
and September 2005 in the H75, H168, and H214 config- 
urations. As summarized in Table [H this combination 
of array configurations provided baseline lengths ranging 
from 30 m to 200 m, and hence sensitivity to structures 
(with uniform brightnesses) as large as about 23'^ Five 
of the six antennas of the ATCA were equipped with 
dual polarization receivers in the 3 mm band, covering a 
frequency range from 85 GHz to 105 GHz. At the wave- 
length of the HCN (1-0) (rest frequency of 88.6318470 
GHz) and HCO+ (1-0) (rest frequency of 89.1885230 
GHz) lines, the primary beam of the antennas is ^36''. 
We configured the correlator to observe the HCN (1- 
0) and HCO+ (1-0) lines simultaneously in two spec- 
tral windows, with the individual windows covered by 
64 channels of width 2 MHz (7 km s~^) each. Under- 
sampling in the lag domain makes the spectral chan- 
nels not completely independent; we therefore binned 
the data into 15 km s~^ channels for making maps. We 
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adopted a pointing center of a2ooo = 23^ 16"^ 11^ and 
^^2000 = — 42°35'05'', coinciding closely with the center 
of the circumnuclear starburst ring. The typical single- 
sideband system temperatures varied over the range 250- 
300 K, and the total on-source integration time was 
about 17 hr. 

We reduced the data in the standard manner using 
MIRIAD, a procedure that involved baseline corrections, 
flagging of bad channels and shadowed baselines, cali- 
bration, and finally imaging. The absolute flux density 
scale was set by Uranus or Mars, excluding baselines for 
which the planets were heavily resolved. We used the 
bright quasar PKS B 192 1-293 for bandpass calibration, 
and the quasar PKS B2326-477 (lying within 5° of our 
target) for complex gain calibration. We also observed 
another quasar, PKS B2255-282 (lying within 15° of our 
target), at frequent intervals to verify the phase solu- 
tions. 

We made channel maps using robust weighting, yield- 
ing a synthesized beam of 2^6 x 2^0, at a PA of 80°. The 
final maps have an rms noise level of 4.5 mJy in HCN (1- 

0) and 4.2 mJy in HCO+ (1-0) for each 15 km s"^ chan- 
nel. Because the circumnuclear emission has a small an- 
gular extent compared with the primary beam, we did 
not apply any primary beam corrections to the maps. 

2.3. Radio observations: Submillimeter Array 

We observed NGC 7552 in ^^CO (2-1) and ^^CO (2-1) 
in the compact configuration of the Submillimeter Ar- 
ray (SMA; Ho etlD l2004D on 4 August 2006. Obser- 
vations were centered at a2000=23^16"^10!7 and (^2000 = 
— 42°35'05V41, coinciding with the center of the circum- 
nuclear starburst ring as defined bv lForbes et al.l (|1994a[ ). 
The primary beam of the antennas at the wavelengths of 
the observed molecular lines is 55'^ Seven out of the 
eight antennas of the SMA were available during our 
observation. As summarized in Table [TJ with baseline 
lengths ranging from 16 m to 69 m, the largest structure 
(of uniform brightness) that we expect to be sensitive to 
is about 16'^ The receiver was tuned and the correlator 
configured to cover the frequency range 227.6-229.6 GHz 
in the upper sideband (USB) thus including the ^^CO (2- 

1) line, and 217.6-219.6 GHz in the lower sideband (LSB) 
thus including the ^^CO (2-1) line. The correlator pro- 
vided a total of 24 spectral windows (chunks) in each 
sideband, with a total of 128 channels and a bandwidth 
of 104 MHz in each chunk. To make the channel maps in 
^^CO (2-1), we averaged 10 channels together resulting 
in a velocity resolution of 10 km s~^. Because the signal 
is weaker in ^^CO (2-1), to make the channel maps in 
this line we averaged 18 channels together resulting in a 
velocity resolution of 18 km s~^. 

We calibrated the data separately for bandpass, ampli- 
tude, and phase using SMA-specific MIR tasks adopted 
from the MMA software package (written in IDL), which 
was originally developed for the Owens Valley Radio Ob- 
servatory (OVRO, Scovill e et al.l il993). The bandpass 
was calibrated using Uranus. For complex gain calibra- 
tion, we observed the quasars J2235-485 and J2258-279 
for a duration of 6 min every 30 min. We used J2235- 
485, which is located at an angular distance of 9? 3 from 
our target source, for phase calibration. J2235-485, how- 
ever, is too weak for amplitude calibration. We therefore 



used J2258-279, which is located at a larger angular dis- 
tance of 15° from our target source, for amplitude cal- 
ibration. Because the amplitude solutions derived from 
the individual scans did not differ significantly between 
adjacent scans, we applied time smoothing scale of three 
scans of J2258-279 to derive more precise amplitude so- 
lutions. We derived the absolute flux calibration scale 
from Uranus, which we assumed had a flux density of 
36.836 Jy at the time of our observation. 

We imaged the data using MIRIAD. To deconvolve the 
sidelobes, we first CLEANed the dirty maps without any 
preselection to judge which features are likely to be real 
and which are likely to be sidelobes. We then placed 
boxes around those features judged to be real on the ba- 
sis that they do not appear to resemble the sidelobes of 
the dirty beam, and that their structure changed con- 
tinuously between adjacent channels. As shown in §3.5[ 
the features in each channel have a sufficiently simple 
structure that we are confident of having picked out only 
those that are real. The ^^CO (2-1) channel maps have 
an angular resolution of T!0 x 2^8 at PA = — 11?9, and a 
rms noise level of 78 mJy beam~^ at a velocity resolution 
of 10 km s~^. The ^^CO (2-1) channel maps were made 
with 18 km s~^ velocity binning, resulting in an angular 
resolution of 6^9 x 2^8 at a PA = —8? 9, and a rms noise 
level of 39 mJy beam~^. Because the features detected 
span an angular extent from the pointing center that is 
significantly smaller than the primary beam, we did not 
apply a primary beam correction to the maps. 

3. RESULTS 
3.1. Optical images 

Figure [1] shows the optical images obtained with Gem- 
ini South. Distinct dust lanes are often seen in late-type 
barred galaxies along the leading edges of the bar, as is 
the case also for NGC 7552. In the g'-band image of Fig- 
ure m we indicate the location of the dust lanes with a 
more direct tracer, dust emission as imaged by Spitzer 
at 5.8 /im. The overall extent of this dust emission is 
indicated by a contour, which perfectly matches the dust 
lanes seen in the optical. 

The Ha image of Figure [T] highlights regions of recent 
activity involving massive stars. The circumnuclear star- 
burst ring as traced in the 3-cm continuum is overlaid on 
the Ha image to show the overall extent of this ring. 
The physical size of the ring at 3 cm is comparable to 
the size of the central bright bulge of NGC 7552 in the 
optical images. A magnified view of the central region 
in Ha and 3-cm continuum is shown in Figure [3l In ad- 
dition to the circumnuclear starburst ring, a bright H II 
region, detected in both Ha and radio continuum, is ap- 
parent ~3.5 kpc east of the nucleus along the bar. This 
H II region was also detected in earlier Ha images by 
iForbes et all (|l994ah . 

3.2. HI 

A map of the integrated H I intensity (zeroth moment) 
is shown in Figure [2] (left panel). As pointed out by 
Dahlem (2005), the H I gas traces the large-scale disk 
of NGC 7552. For ease of comparison, the optical i'- 
band image of the galaxy shown earlier in Figure [T] is 
overlaid as contours. Along the major axis the H I disk 
can be detected out to a projected radius of about 2100'^ 
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(~20 kpc), about twice as far out as the optical disk. 
Similarly, along the minor axis, the H I disk can be de- 
tected out to a projected radius of about 121'^ (~13 kpc), 
again about twice as far out as the optical disk. There 
is a prominent ring-like enhancement in H I that peaks 
at or close to the ends of the large-scale optical bar. Fi- 
nally, a central depression can be seen in the integrated 
H I intensity map. This central depression has a size 
that closely coincides with the angular extent of the cir- 
cumnuclear molecul ar ga s detected in our observations 
as described below ( §3.41 and §3.5p . 

Figure [2] (right panel) shows a map of the intensity- 
weighted mean H I velocity (first moment). Assuming 
a systemic velocity of 1586 km s~^ as found in the opti- 
cal and as derived below (^33j) from our molecular line 
observations, emission can be seen spanning the velocity 
range from —100 km s~^ to 100 km s~^ about the sys- 
temic velocity. As can be seen, the kinematic major axis 
is aligned roughly east- west, coincident with the major 
axis of the optical bar. Blueshifted emission (negative 
velocities) originates primarily westwards and redshifted 
emission (positive velocities) primarily eastwards of the 
galactic center. On close scrutiny, the position angle of 
the kinematic major axis can be seen to change with 
galactic radius, creating an S-shaped pattern in the gas 
kinematics. Such S-shaped kinematic patterns are com- 
monly attributed to non-circular motions under the non- 
axisymmetric gravitational potential of a galactic bar 
(see §5.1l for detailed discussion). Far beyond the optical 
bar however, the position angle of the kinematic major 
axis continues to change, suggesting that the outskirts of 
the galaxy may be warped (by gravitational interactions 
with one or more neighboring galaxies as described in 

m- 

3.3. Radio Continuum 

Figure [3] shows the central part of our radio contin- 
uum images at 6 cm (left panel) and 3 cm (right panel), 
both of which were recreated from the data taken by 
iForbes et al.l (jl994al ). Over the field shown, the radio 
continuum emission originates from a ring, which as men- 
tioned by Forbes et al.i (; 1994 a) corresponds to the cir- 
cumnuclear starburst ring, as can be seen in Figure [3] 
(right panel) where the 3-cm image (contours) is over- 
laid on our Ha image (color). This ring has an outer ex- 
tent comparable to the central H I depression mentioned 
above. 

Five peaks are clearly visible in the 3-cm image (right 
panel), whereas only three peaks are visible in the 6-cm 
image (left panel) that has a lower angular resolution 
than the 3-cm image. To study the spectral index of 
each knot, we have convolved the 3-cm image to the same 
angular resolution as the 6-cm image. After degrading 
the angular resolution, we find that only three knots are 
discernible in the-3 cm image, which now resembles the 
6-cm image. The spectral indices (a, where S (x u'^ with 
S the fiux density and u the observing frequency) of the 
three knots (see the left panel of Fig. [3]) are —1.11 ±0.02 
(upper right), -0.84 ±0.03 (upper left), and -0.72 ±0.05 
(bottom left), with an averag e value of —0.89 ±0.0 6 that 
is similar to that inferred by iForbes et al.l (|1994a[ ) . The 
negative spectral indices indicate that the emission of all 
the knots is dominated by, if not produced entirely by, 
(non-thermal) synchrotron emission. 



3.4. HCNandHCO"' 

Figure 2] shows our channel maps of the HCN (black 
contours) and HCO+ (1-0) (cyan contours) lines. Both 
HCN (1 - 0) and HCO+ (1 - 0) emission are concen- 
trated in two knots located to the north and south of the 
galactic nucleus. The northern knot is detectable over a 
velocity range from —133 km s~^ to 32 km s~^ with re- 
spect to the optically-determined systemic velocity of the 
galaxy, and is therefore preferentially blueshifted. The 
southern knot is detectable over a velocity range from 
about —43 km s~^ to 92 km s~^ with respect to the sys- 
temic velocity, and is therefore preferentially redshifted. 
The centroids of both knots shift from west to east with 
increasing (i.e., from blueshifted to redshifted) velocities. 

Figure [5] (left panel) shows a map of the integrated 
HCN (1-0) intensity in black contours. The overall 
size of the emitting region is comparable with the cen- 
tral H I depression described above in §3.2[ suggesting 
that the gas density is so high here that the atomic hy- 
drogen gas has been largely converted to molecular hy- 
drogen gas. The HCN (1-0) contours is overlaid on a 
greyscale image of the 3-cm continuum emission shown 
earlier in Figure [3] (right panel) but now convolved to 
the same angular resolution as the HCN (1-0) map. At 
this angular resolution, the 3-cm continuum map shows 
two knots just like the HCN (1 - 0) map. The observed 
HCN (1-0) emission therefore originates from the cir- 
cumnuclear starburst ring, and exhibits a central depres- 
sion or possibly a hole. Notice the small displacement 
between the centroids of the HCN (1-0) gas and radio 
continuum knots, a point we shall return to in §6.21 

Figure [5] (right panel) shows a map of the intensity- 
weighted mean HCN (1-0) velocity. A velocity gradient 
is discernible roughly along the east-west direction, with 
the kinematic major axis lying at a position angle of 110°. 
The kinematic major axis in HCN (1 - 0), as well as 
in HCO^ (1 - 0) as described below, are both aligned 
with the kinematic major axis at the innermost region of 
the large-scale H I disk as described in §3.21 The good 
agreement suggests that there is a continuous velocity 
field between the outermost region in molecular gas and 
innermost region in neutral atomic gas. 

Figure [6] shows our moment maps of the HCO^ (1 - 

0) line. The emission in all these maps is virtually iden- 
tical to that seen in HCN (1-0). Both HCN (1 - 0) 
and HCO+ (1-0) have comparable critical densities of 
~10^ cm~^, suggesting that both are tracing the same 
regions of dense gas. 

3.5. CO 

Figure [71 shows our channel maps of the ^^CO (2 - 

1) line, which by contrast with both the HCN (1-0) 
and HCO+ (1-0) has a much lower critical density of 
just '^10^ cm~^. Only a single compact component is 
discernible that shifts from west to east across the cen- 
ter of the galaxy with increasing velocity, just like that 
seen in HCN (1 - 0) and HCO+ (1 - 0) as described 
above in §331 We have convolved both the HCN(l-O) 
and HCO+(1-0) maps to the same angular resolution 
as the ^^CO (2-1) map, and found that all show the 
same spatial- kinematic structure suggesting that they all 
trace the same global features. In the ^^CO (2-1) map, 
emission is detectable over the velocity range from about 
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—95 km s ^ to 100 km s ^ with respect to the systemic 
velocity, comparable to the detectable velocity range of 
the emission in both the HCN (1 - 0) and HCO+ (1 - 0) 
maps. 

Figure [8] (left panel) shows the integrated ^^CO (2 - 1) 
intensity map in black contours overlaid on the integrated 
HCN (1-0) intensity map, taken from Figure [5] but con- 
volved to the same angular resolution as the ^^CO (2 - 
1) map, in grayscale. Figure [8] (right panel) shows a map 
of the intensity- weighted mean ^^CO (2-1) velocity. As 
can be seen, the compact central component detected in 
-'^^CO (2-1) has a similar size and kinematics as the cir- 
cumnuclear HCN (1 - 0) (and HCO+ (1 - 0)) emission. 
In Figure m we show position- velocity (P-V) diagrams of 
the HCN (1 - 0), HCO+ (1 - 0), and ^^CO (2 - 1) emis- 
sion along a position angle of 110°, which corresponds to 
the orientation of the kinematic major axis. As can be 
seen, the P-V diagram in ^^CO (2-1) closely resembles 
that in the higher density tracers of molecular hydrogen 
gas. We therefore believe that the observed ^^CO (2-1) 
emission originates from essentially the same regions of 
the circumnuclear starburst ring as the observed HCN (1 

- 0) or HCO+ (1-0) emission. Of course, any central 
depression or hole of the same size as that seen in HCN (1 

- 0) and HCO+ (1-0) would not be discernible at the 
angular resolution of our ^^CO (2-1) observation. 

The P-V diagrams of Figure [9] reveal that the systemic 
velocity at the center of the circumnuclear starburst ring, 
which presumably coincides with the dynamical center of 
the galaxy, is about 1580 km s~^. This is consistent with 
the optically-determined value of 1586±5 km s~^ (RC3), 
but significantly smaller than that determined from the 
global H I p rofile of 1608 ±5 km s"^ (The H I Parkes Ah 
Sky Survev: [Barnes et al.l[200Th . 

Figure [To] shows our channel maps of the ^^CO (2-1) 
line. By contrast with the simplicity of the ^^CO (2 - 
1) channel maps, three components are now discernible: 
a compact central component that is the counterpart of 
that seen in ^^CO (2 - 1), and two extensions on oppo- 
site sides of the center with one arm extending east from 
north and the other arm west from south. The com- 
pact central component is detectable over the velocity 
range from about —125 km s~^ to 135 km s~^ with re- 
spect to the systemic velocity, somewhat larger than the 
detectable velocity range of its counterparts in ^^CO (2 

- 1), HCN (1 - 0), and HCO+ (1 - 0). By contrast, 
both the eastern and western arm-like components are 
detectable over a narrower velocity range than the com- 
pact central component in any of the observed molec- 
ular lines. The western arm is detectable over the ve- 
locity range from about —85 km s~^ to 55 km s~^ and 
is therefore preferentially blueshifted, whereas the east- 
ern arm is detectable over the velocity range from about 
—55 km s~^ to 75 km s~^ and is therefore preferentially 
redshifted. 

Figure [TT] (upper left panel) shows a map of the in- 
tegrated ^^CO (2-1) intensity in black/white contours 
overlaid on the optical z'-band image of the inner part of 
the galaxy. Dust lanes are visible as silhouettes in the 
optical image running east and west of the center, as well 
as in emission at 5.8 fim as s hown by the cyan con tour 
based on the observations by iKennicutt et al.l ()2003i) in 
their Spitzer SINGS project. The eastern and western 



^^CO (2-1) arms coincide with the innermost regions 
of these dust lanes. These arms appear to connect with 
the circumnuclear starburst ring at or near the locations 
where we detected two knots of dense molecula r ga s as 
traced in HCN (1 - 0) and HCO+ (1 - 0) (see g331), as 
shown in the lower left panel of Figure [TT] where the inte- 
grated ^^CO (2-1) intensity is plotted in green contours 
and integrated HCN (1-0) intensity in black contours. 
We note that, for the same ratio in intensity between the 
arms and compact central component, the stronger west- 
ern arm should be detectable (at 4.5cr) at the sensitivity 
of our HCN (1-0)/HCO+ (1-0) maps. 

Figure [TT] (lower right panel) shows a map of the 
intensity- weighted mean ^^CO (2-1) velocity. At our 
angular resolution, the kinematics of the compact cen- 
tral component cannot be completely separated from the 
individual kinematics of the eastern and western arms, 
and vice versa. Nevertheless, the S-shape velocity field 
(contours) indicates the existence of non-circular motion, 
which is similar to what is observed in other barred galax- 
ies. In Figure [9] we show a P-V diagram made from the 
CO (2-1) channel maps at a position angle of 110°. 
As can be seen, the kinematics of the compact central 
component are similar to its counterparts in the other 
observed molecular lines. The western arm is preferen- 
tially blueshifted whereas the eastern arm preferentially 
redshifted, and therefore both arms share the same over- 
all kinematics as the large-scale galactic disk seen in H I. 

We derive the mass of molecular hydrogen gas in all 
the three features detected in ^^CO (2-1) in the stan- 
dard manner using the Galactic conversion factor be- 
tween CO intensity and mass of molecula r hydrogen gas 
ofXc o = 3 X 10^0 cm (K km s'^)'^ (jSolomon et al.l 
Il987[ ). This conversion is based on the ^^CO (1 - 0) 
line, and in the present situation requires knowledge 
of th e ratio in ^^CO (2 - 1) to ^^CO (1 - 0) intensi- 
ties. lAalto et al.l ([T995I ) find a ratio in brightness tem- 
peratures between these two lines of 1.2 ± 0.1 based 
on single-dish observations with the Swedish-ESO 15- 
m Submillimeter Telescope (SEST) towards the center 
of NGC 7552. For simplicity, we shall therefore assume 
a ratio in brightness temperature of unity for ^^CO (2 

- 1) to ^^CO (1 - 0) throughout this work. Accord- 
ingly, the compact central component (corresponding to 
the circumnuclear starburst ring, and perhaps regions 
within), which has an integrated intensity in ^^CO (2 

- 1) of 1487.1 ± 20.3 Jy km s~^, has a gas mass of 
(3.2±0.1)xl0^ Mo, including the mass of the heavy ele- 
ments other than molecular hydrogen (Mgas = I.SGMyl^)- 
The integrated ^^CO (2-1) intensity detected over the 
entire region is 2330.5 ± 34.7 Jy km s~^, and hence the 
mass of molecular gas in the eastern and western arms 
combined is (1.7 ± 0.2) x 10^ M©. 

The dense gas fraction can be estimated by compar- 
ing the mass derived from ^^CO (2-1) and HCN (1 

- 0). Assuming that the gas traced by HCN (1 - 
0) is gravitationally bound (jSolomon et al.l ri992f ). the 
HCN-to-H2 conversion factor Xhcn = ^(H2)/I/'hcn 
is -20 Mq (K km s"^ pc~^)~\ where A^(H2) is the 
column density of molecular hydrogen gas and L'hcn 
the HCN (1-0) luminosity. For the measured L'hcn ^ 
7.5 X 10^ K km s~^ pc~^, the corresponding mass in 
dense gas, as traced in HCN (1-0), is -1.5x10^ M©. The 
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dense gas fraction is therefore about 50% at the compact 
central component, much higher than is usually found 
in Galactic giant molecular clouds. We note, however, 
that the inferred dense gas fraction may contain consid- 
erable uncertainty because the conversion factors Xhcn 
and Xco used may not be appropriate for the circum- 
nuclear regions of galaxies (although the use of different 
conversion factors can either raise or lower the dense gas 
fraction). For example, based on virial arguments for the 
circumnuclear starburst ring of the nearby barred spiral 
galaxy IC 342, Meier & Turner (2001) argue that Xco 
is only 1 x 10^^ cm ~^ (K km s~^)~^ or even smaller, in 
which case the dense gas fraction at the compact central 
component becomes nearly 100%. 

4. PHYSICAL PROPERTIES OF THE 
CIRCUMNUCLEAR MOLECULAR GAS 

In this section, we derive the physical properties of the 
molecular gas in the circumnuclear starburst ring based 
on the ratio in brightness temperatures of the molecu- 
lar lines observed. We then investigate the dynamical 
stability of this gas. 

4.1. Emission Recovered in our Interferometric 
Observations 

Line ratios inferred with an interferometer can be com- 
promised if there exists uniformly-bright features with 
sizes beyond tha t detec table even by the shortest base- 
line. lAalto etld] ()1995l ) have observed the central region 
of NGC 7552 with the SEST in ^^CO (1 - 0), ^^CO (2 - 
1), i^CO (1 - 0), ^^CO (2 - 1), and HCN (1 - 0). They do 
not show line profiles for any of the transitions observed, 
but only quote the line intensities or ratio in brightness 
temperatures of the different lines. After convolving our 
channel maps to the same angular resolution as in the 
observation of the corresponding line with the SEST, we 
found that our SMA observation recovered 54% ± 9% of 
the ^^CO (2 - 1) emission detected by the SEST. Given 
typical systematic uncertainties of about 20% in the ab- 
solute flux calibration at millimeter wavelengths, we may 
well have recovered the bulk of the circumnuclear emis- 
sion in ^^CO (2 - 1). Our ATCA observation recovered 
76% ± 16% of the HCN (1 - 0) emission detected by 
the SEST, suggesting once again that the bulk, if not 
all, of the circumnuclear emission in this line has been 
recovered. Surprisingly, our SMA observation recovered 
only 26% ± 10% of the ^^CO (2 - 1) emission detected 
by the SEST, despite the fact that ^^CO (2 - 1) traces 
higher column densities and therefore more compact re- 
gions than ^^CO (2 - 1) as evidenced by the maps pre- 
sented in Figure [HI The ^^CO (2 - 1) line however was 
detected at a confidence level of only 3-4cr in the SEST 
observation, and we therefore do not place a high degree 
of confidence on this particular comparison. 

4.2. Physical Properties 

We computed the physical properties (density and tem- 
perature) of the molecular gas associated with the cir- 
cumnuclear starburst ring from the ratio in brightness 
temperatures of the HCN (1 - 0) to ^^CO (2 - 1) lines 
using the Large Velocity Gradient (LVG) approxima- 
tion (e.g., iGoldreich fc KwanI [1971 [Scoville Solomod 
[1971 . We used the ^^CO (2 - 1) line because it IS more 



likely than the ^^CO (2-1) line to trace the same region 
as the HCN (1 - 0) line, as the maps we presented earlier 
indicate (see §3.51 and Figure [8]). HCN (1 - 0) is a tracer 
of high density molecular hydrogen gas; as the ^^CO (2 - 
1) emission is much more optically thin than the ^^CO (2 
- 1) emission, the former is more strongly weighted to- 
wards regions of high column density and therefore also 
density. Indeed, as we shall show, the -'^^CO (2-1) emis- 
sion is optically thin and therefore a good tracer of col- 
umn density. 

To derive the ratio between the HCN (1-0) and 
^^CO (2-1) intensities, we first convolved our HCN (1 - 

0) map to the same angular resolution as our ^^CO (2 - 

1) map (see Figure [8]). The circumnuclear starburst ring 
is spatially unresolved at this angular resolution, and so 
the line ratio we derive for HCN (1 - 0) to ^^CO (2 - 
1) of 2.0 ± 0.2 is an average over the entire ring. We 
assumed Galactic abundances for the ^^CO and HCN 
molecules of [^^C0]/[H2] ^ 10"^ (jSolomon et al.l [T979I ) 
and [HCN]/[H2] ^ 2 x 10"^ (jlrvine et al.[[T987l ). Finally, 
we estimated a velocity gradient of 0.24 km s~^ pc~^ 
from the P-V diagram along the kinematic major axis 
of the HCN (1-0) emission as shown earlier in Figure [H 

Figure [12] shows a plot of the density of molecu- 
lar hydrogen gas, n(H2), as a function of the product 
Z{^^CO)/{dv/dr), where Z{^^CO) is the ^^CO abun- 
dance and dv/dr the velocity gradient along the radial 
direction. For the adopted or inferred values above, 
Z{^^ CO) /{dv/dr) ^ 4 x 10-^ which is displayed with 
the blue vertical line in each panel. The different pan- 
els in this figure correspond to different kinetic tem- 
peratures (from 10 K to 1000 K). The red contours in 
each panel indicate acceptable solutions in the n(H2) and 
Z(^^ CO) /{dv/dr) phase space for different HCN (1 - 0) 
to ^^CO (2-1) line ratios based on our LVG calcula- 
tions. The black lines indicate the ^^CO (2-1) opacity. 
As can be seen, for a HCN (1 - 0) to ^^CO (2 - 1) line 
ratio of 2.0 ± 0.2, and Z{^^ CO) /{dv/dr) ^ 4 x 10-^ so- 
lutions can only be found for kinetic temperatures well 
above 100 K, for which the ^^CO (2 - 1) emission is op- 
tically thin. At such temperatures, the density of the 
molecular gas (averaged over the line-emitting regions) 
is around 10^^^ cm~^. A similar result has been seen 
in the circumnucle ar starburst ring of the barred spiral 
galaxy NGC 1097 (iHsieh et al.ll2008[ ). 

On the other hand, the adopted Galactic abundance for 
i^CO of [^^C0] /[H2] ^ 10-^ may be too high for star- 
burst regions (faniguchi fc Ohvamal [1998^ . The ^^CO 
abundance is controlled primarily by fractionation reac- 
tions and photodissociation. Fractionation reactions en- 
hance the ^^CO abundance, but such reactions are only 
effective in regions where t he temperature is less than 
~35 K ([Watson et al.[ [19761 ) and hence not exposed to 
vigorous star formation. Photodissociation, naturally, is 
even more effective in regions where the star-formation 
rate is higher. If the combination of these factors reduces 
[^^C0]/[H2] by, say, an order of magnitude to ~10~^, 
then acceptable solutions to the observed HCN (1-0) 
to ^^CO (2-1) line ratios can be found for temperatures 
^50-100 K and densities ~10^ cm"^. Thus, unless the 
^^CO abundance is even lower than the range of values 
considered above, the molecular gas that we detected 
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in the starburst ring is relatively warm and dense, as 
would be expected for gas associated with regions that 
have recently undergon e vigo rous star formation (see also 
[Matsushita et al. 19981, |2Q10D. 

In many galaxies exhibiting nuclear activity, it is often 
difficult tell whether this activity is produced by activ e 
galactic nuclei (AGN) or starburst. iKohno et al I ()2QQ1I) 
suggest that such galaxies can be distinguished by plot- 
ting the line ratio of HCN (1 - 0) to HCO+ (1 - 0) 
against that of HCN (1 - 0) to ^^CO (1 - 0) (Figure ??). 
They show that galaxies dominated by or exhibiting only 
AGN activity lie in the upper right corner of such a di- 
agram, having high values in both line ratios. By com- 
parison, galaxies exhibiting only starburst activity or a 
mix of starburst and AGN activity lie in the lower left 
corner, having low values in both line ratios. In the case 
of NGC 7552, we can clearly see that the nuclear activ- 
ity is dominated by a starburst (i.e., circumnuclear star- 
burst ring), and our observation are able to confidently 
separate the molecular gas associated with the starburst 
from any associated with the nucleus in HCN (1-0) and 
HCO+ (1 - 0), and also to some degree of confidence in 
i^CO (2 - 1) and ^^CO (2 - 1) (based on the similarity 
of the gas kinematics in all four lines). Such is not the 
case fo r the galaxies used to make the abovementioned 
plot in IKohno et all (2001). For NGC 7552, the HCN (1 

- 0) to HCO^ (1-0) line ratio, averaged over the cir- 
cumnuclear starburst ring, is 1.13 ± 0.04. The HCN (1 

- 0) to ^^CO (1 - 0) ratio, assuming that the ^^CO (1 - 
0) is equal to the ^^C O (2 - 1) intensity for the reason 
mentioned earlier (see §4.ip . is 0.181 ±0.005. This places 
NGC 7552 firmly in the region occupied by galaxies with 
nuclear activity dominated b y a starburst, enhan cing the 
validity of the plot made by iKohno et al.l ()2001f ) . 

4.3. Dynamical Stability 

We examine whether the molecular gas in the circum- 
nuclear starburst rin g is dynamical ly stable according to 
the Toomre criterion ( Toomrelll964f ). For a disk to be sta- 
ble against collapse due to its own gravity, the Toomre 
Q parameter, which is defined as the ratio of the criti- 
cal surface gas density to the actual surface gas density, 
Scrit/Sgas, is required to have a value equal to or exceed- 
ing unity. 

We compute the average surface gas density in a cir- 
cular region of radius of 10^' (~1 kpc), just encompass- 
ing the circumnuclear starburst ring, using the molecular 
gas mass derived in §3.51 to find Sg^s ^ 900 Mo pc-^ 
Because the ^^CO (2 - 1) emission likely exhibits a cen- 
tral depression or hole within the circumnuclear starburst 
ring, the value we computed is likely to be a lower limit. 
We estimated th e critical surf ace gas density using the 
relationship from IKohno et all (jl99£) 

= (l5l^) ( o.57km%pc-0 

(1) 

where a is the one-dimensional velocity dispersion, for 
which we estimated a value of 12 km s~^ based on 
a second moment map of the ^^CO (2-1) emission 
(not show n), a nd the epicyclic frequency. As shown 
below in §5.2| we derived an epicyclic frequency of 
0.4 km s~^ pc~^ at the radius of the circumnuclear star- 



burst, and therefore Scrit 

^ 354 Mo pc-^ The Toomre 
Q parameter is therefore no larger than 0.4 in the cir- 
cumnuclear starburst ring, indicating that the molecular 
gas there is (highly) unstable to gravitational collapse. 

5. GLOBAL KINEMATICS 

In this section, we derive the global kinematics (ro- 
tation curve) of NGC 7552 and, assuming a pattern 
speed for the galactic bar, estimate the predicted lo- 
cations of dynamical resonances. We then compare 
these predictions with the locations of prominent fea- 
tures in NGC 7552 that in other disk galaxies are com- 
monly attributed (but seldom demonstrated explicitly 
to correspo nd) to dynani i cal reso nances; in the case of 
NGC 7552, iForbes et al.l (jl994br ) have highlighted fea- 
tures that possibly correspond to dynamical resonances. 
Our overarching goal is to determine whether we can un- 
derstand the global kinematics of NGC 7552 sufficiently 
well so as to test ideas for the formation of dense molec- 
ular gas and stars in its circumnuclear starburst ring as 
discussed in the next section. 

5.1. Rotation Curve 

As pointed out in ^3.21 the intensity- weighted mean H I 
velocity map (Fig. [2]) exhibits an S-shaped asymmetry 
characteristic of periodic orbits predicted theoretically 
for barred galaxies (e.g.. iAthanassoulalll992l and refer- 
ences therein). In general, such orbits are elliptical (and 
centered on the dynamical center of the galaxy) , and have 
their major axes either parallel {xi orbits) or orthogonal 
{x2 orbits) to the major axis of the bar, changing from 
one to the other between the locations of dynamical res- 
onances. Specifically, beyond the ILR (inner Lindblad 
resonance) out to corotation (CR), orbits are xi. Within 
the ILR, orbits are X2. If there is both an outer ILR 
(oILR) as well as an inner ILR (ilLR), orbits are xi inside 
the ilLR and X2 between the ilLR and oILR. Between CR 
and the outer Lindblad resonance (OLR), orbits are xi. 
Close to and beyond the major axis of the bar (predicted 
theoretically to terminate at or before CR), both xi and 
X2 orbits become progressively more (and are predicted 
to be usually close to) circular. Close to a bar potential, 
xi is the predominant orbit, hence the major kinematic 
axis is closely aligned with the galactic bar. The change 
from xi to X2 orbits and back between dynamical reso- 
nances causes the S-shaped twist in the major kinematic 
axis on scales up to that comparable with the length of 
the galactic bar. 

To accurately derive the rotation curve of NGC 7552 
from our measurements of the projected galaxy kinemat- 
ics, we require a priori knowledge of the axial ratios 
(which change with galactic radius) of the xi and X2 
orbits, and the inclination of the galaxy (including any 
variation in inclination with galactic radius if the disk is 
warped). As these parameters (apart from perhaps incli- 
nation) are not in general known, we assumed a constant 
inclination of 25° for NGC 7552 based on the axial ratio 
of its stellar disk as described in §1. For the major kine- 
matic axis of the disk, we used a position angle of 110° as 
indicated by the first moment map in all four molecular 
lines as well as the inner region of the first moment map 
in the H I line. 

The HCN (1 - 0) to HCO+ (1 - 0) emissions orig- 
inate primarily, if not entirely, from the circumnuclear 
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starburst ring, and are therefore useful for inferring the 
rotational velocity of this ring. Sampling only a limited 
and furthermore spatially unresolved radial range, how- 
ever, these lines do not provide a meaningful rotation 
curve at the center of the galaxy. The ^^CO (2 - 1) 
line samples a larger radial range from (presumably) the 
ring to the innermost regions of the dust lanes, and is 
useful for making a crude estimate of the rotation curve 
over this range. The H I line provides a relatively robust 
estimate of the rotation curve beyond the outer regions 
probed in ^^CO (2 - 1) to far beyond the stellar optical 
disk. 

The rotation curve is directly derived from P-V dia- 
grams of ^^CO (2-1) and H I from cuts along 110°, 
as shown earlier in Figure [9l with the enve lope-tracing 
meth od (see details and more examples in iSofud 1 1 9961 
11997). The method traces terminal velocity (or tangent- 
points) on PV diagram and corrects for the ISM velocity 
dispersion. The deprojected rotational velocity is derived 
from the terminal velocity Vt'. 



Vrot =Vt/smi - (cr^bs 



(2) 



where i is the inclination, dobs the velocity resolution of 
the observation, and (Jism the velocity dispersion of in- 
terstellar gas. We adopt an ISM velocity of 7 km s~^ 
fetark fc Brand 1989: Malhotra 1994). The terminal ve- 
locity is the velocity at which the intensity {It) is equal 
to: 

It= \{vlm..f + (3) 



on the P-V diagram, where /max and lie are, respec- 
tively, the maximum intensity and intensity of the low- 
est (typically the 3a) contour level. If the intensity is 
high enough. It ~ I\c and is therefore approximately de- 
fined by the locus of the rotation curve, r] is taken as 
0.2 as suggested in "Sofue (1996), corresponding to the 
20% level of the intensity profile at a fixed position, i.e., 
It c± 0.2 X /max- The accuracy of the obtained Vt along 
with the Vrot by the envelope-tracing method is ~ ±15 
km s~^, correspondin jS!: to y ±15 km s~^/sinz in depro- 
jected velocity (SoSi[l996l [1993). 

In Figure [T4j we plot the projected rotation curve de- 
rived from the ^^CO (2 - 1) (from to - 0.8 kpc) and H I 
(beyond 1 kpc) lines onto the P-V diagram of both these 
lines cut along a position angle of 110°. The projected 
uncertainties of ±15 km s~^ are plotted on the approach- 
ing side of the rotation curve to indicate the size of the 
error bars. Due to the imperfect physical connection be- 
tween the two tracers and asymmetry of P-V diagrams, 
the rotation curve between 0.5 to 1 kpc is poorly con- 
nected; we indicate these regions with grey shading in 
Figure [TH 

The deprojected inner and the entire derived rotation 
curves are shown respectively in Figures ?? and ??. The 
open and solid black circles denote the receding and ap- 
proaching sides of the rotation curve respectively, and 
the red open circles indicate the average of the two sides 
at the radii where both receding and approaching veloc- 
ities are seen. The deprojected uncertainties are plotted 
in the lower right of each panel. 

5.2. Dynamical Resonances 



Because the intensity of the approaching side in the 
H I map is significantly weaker than that of the reced- 
ing side, together with the asymmetry of the derived ro- 
tation curve on each side, we determine the dynamical 
resonances for each side of the galaxy separately. We 
compute where dynamical resonances are theoretically 
predicted to occur using the usual assumption that the 
angular velocity of the bar (its pattern speed, Qp) is equal 
to the angular velocity of the galaxy at the outermost 
ends of the bar (i.e., the CR). We determined the length 
of the bar in NGC 7552 from the Two Micron All-Sky 
Survey (2MASS) image at 2.2 /im, assuming that the bar 
ends where the two prominent spiral arms first appear to 
emerge at a radius of about 5 kpc. 

The average rotation curve shows that the angular ve- 
locity of the galaxy at the bar's end (5 kpc) is about 
38 km s~^ kpc~^, which we therefore take to be the pat- 
tern speed of the bar. In Figure [TH we plot the angular 
velocity (velocity divided by radius, V/R) of the galaxy, 
Q (filled circles), as derived from the rotation curve. The 
soHd line is an interpolation between adjacent angular ve- 
locity data points, from which we derived the epicyclic 

frequency, a^, of the orbits, where tz'^ = 2 ± 

We then computed the Vt — n/2 curve (triangles) to infer 
the location(s) of the ILRs, the Vt — n/A curve (diamonds) 
to infer the location of the Ultra Harmonic Resonance 
(UHR), and the Vt^ n/2 curve (squares) to infer the lo- 
cation of the OLR. The uncertainties in these curves are 
rather large so we do not show the error bars in Fig- 
ure [16] for clarity, but show the possible radial ranges of 
the resonances taking these uncertainties into account. 

We therefore constrain the following possible locations 
of dynamical resonances in NGC 7552. The ^ n/2 
curve (open squares) cross the value of Vt^ at >8 kpc, 
the predicted location of the OLR. Theory predicts that 
the OLR in spiral galaxies should occur just beyond the 
outermost reaches of their spiral arms, which in the case 
of NGC 7552 can be traced out to a (deprojected) radius 
of about 10 kpc. The computed Vt — n/A curve implies 
that the UHR occurs at roughly 2-4 kpc. Forbes et al.l 
(1994b) suggest that an isolated bright H II region lo- 
cated at a (deprojected) radius of 3.5 kpc coincides with 
the UHR. The Vt - k/2 curve implies two ILRs. The 
oILR is located at a radius of about 1.7 kpc, which is be- 
yond the circumnuclear starburst ring that has a radius 
of about 0.5 kpc. Further inwards the angular velocity 
might eventually drop to a value of km s~^ (in prin- 
ciple) at the center of the galaxy, and so an ilLR could 
exist somewhere within the circumnuclear starburst ring. 
On the other hand, we do not exclude the possibility that 
the rotation curve keeps rising due to a black hole in the 
galactic center, resulting in the absence of an ilLR. Our 
finding that the circumnuclear starburst ring is located 
within the oILR and beyond, pos sibly, the ilLR is consis - 
tent with theoretical predictions (|Buta fc Combeslll996f ). 
The reasonable agreement between the computed and an- 
ticipated locations of dynamical resonances in NGC 7552 
provides a measure of confidence in our understanding of 
the global kinematics of this galaxy. 



6. STAR FORMATION IN THE CIRCUMNUCLEAR 

RING 
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In this section, we investigate the star formation prop- 
erties of the circumnuclear starburst ring, determine 
whether the observed propert ies of the ring are co nsistent 
with the model proposed bv lKohno et al l ()1999f ). and if 
not what modifications are required. 

6.1. Star Formation Rate 

The radio continuum observations of iForbes et al.l 
(|1994al ), the data from which we reprocessed here, pro- 
vides the cleanest (i.e., not subject to extinction, emis- 
sion detected at high S/N) measure of the star formation 
rate (SFR) in the circumnuclear ring. Like Forbes et al. 
([l994a), we us ed the standard prescription suggested by 
iCondonI (|1992f ) to derive the supernova rate as ~0.1 yr~^ 
from the measured radio luminosity for the ring. We then 
estimated the SFR fr om the su pernova rate in the man- 
ner also suggested bv ICondonI (|19 92). but corrected for 
a Salpeter Initial Mass Function (IMF) as suggested by 
ICondon et all (|2QQ2D . as about 15 M© yr^^ An inde- 
pendent estimate of the SFR in the ring can be made 
from hydrogen recombination lines in the optical and 
near-infrared, although such estimates are affected by 
extinction. The luminosity of the Ha line of 1.5 x 10^^ 
erg s~^ derived from the Gemini im age, using the re - 
lationship of Lua - SFR suggested by ^Ke nnicuttl (| 19981), 
i mplie s a SFR of about 12 M© yr~^. Mo or wood fc Olival 
(|1988D find a luminosity in the Br7 line of 9.8 x 10^^ erg 
s~^ for the circumnuclear ring. Using the relationship 
between the Br-7 line luminosity and SFR suggested by 
iKennicuttI (|1998f ). the SFR is about 8 M© yr-^ The cor- 
relation between the global HON (1 - 0) and far- infrared 
luminosities found bv^ Gao & Solomonl (2004) for normal 
to extreme star-forming galaxies can be used to provide 
an estimate of the SFR using the measured HCN (1-0) 
luminosity (i^^^j^) as a proxy for the infrared luminosity 
(Lir). With an L'^q^ ^ 7.5 x 10^ K km s"^ pc^ and 
hence Ljr of (5.5 ±0.1) x 10^^ in the circumnuclear 
starburst ring, the corresponding star formation rate (for 
a Salpeter IMF) is 11.0 ± 0.2 M© yr-^ Therefore, in 
summary, the SFR in the circumnuclear starburst ring 
of NGC 7552 is about 10-15 M© yr-^. 

For comparison, the global SFR derived in the same 
manner from observations of the entire galaxy in the ra- 
dio continuum at 20 cm bv lSchmitt et al.l (|2006f ) is about 
18 M0 yr~^. The global star formation rate derived from 
the Ha line, based on our Gemini observations, is about 
47 M0 yr~^. The global star formation rate can also be 
estimated fro m the far-IR ni e asurem ents of NGC 7552 
with IRAS by ISanders et all (|2003l) . They find an in- 
frared luminosity for this galaxy of 1.1 x 10^^ L©, thus 
placing NGC 7552 right at the threshold for a Lu- 
minous Infrared Galaxy. Using as above the conver- 
sion between inf r ared l uminosity and SFR suggested by 
IGao fc Soloi^ (f200l for a Salpeter IMF, the SFR is 
about 22 M0 yr~^. Because the SFR estimated from Ha 
can be highly uncertain (e.g., up to 99% of UV photons 
can leak out of H II regions; iKennicu 1[^98), we give 
more weight to the other methods of estimating SFR 
as described above. Thus, the circumnuclear region of 
NGC 7552 contributes ~50 - 80% of the galaxy's entire 
star formation, and being strongly concentrated in the 
circumnuclear ring, results in an extremely high star for- 
mation rate surface density of >4 x 10~^ M© yr~^ pc~^. 



The SFR in the circumnuclear ring, divided by the 
mass of molecular gas present, is about (3-5) x 10~^ 
yr~^. The inverse of this ratio represents the timescale 
for gas exhaustion and is 3 x 10^ yr. Such a short gas 
exhaustion timescale is comparable to that of interacting 
objects and ULIRGs (few x 10^ yr) and about 10 time s 
shorter than that in normal galaxies (|Bigiel et al.ll2008[ ). 
Note that even though the gas exhaustion timescale of 
circumnuclear starburst ring is short, as the galaxies with 
circumnuclear rings are generally harboring stronger bars 
than the galaxies lacking rings, the gas is likely to be 
nearly continuously replenished at the circumnuclear re- 
gion through, presumably, the bar-driven gas infiow. 

6.2. Formation of Dense Molecular Gas and Stars at 
the Circumnuclear Ring 

The most prominent optical features — circumnuclear 
starburst ring, large-scale bar, and spiral arms — in 
NGC 7552 are remarkably similar to those of NGC 6951, 
including even the orientation of these features on the 
sky, thus providing a simple and direct comparison. Like 
NGC 7552, the bar in NGC 6951 is oriented east-west, 
with spiral arms emerging from the bar and curling north 
on the western side and south on the eastern side. Dust 
lanes that run along the bar intersect the circumnuclear 
starburst ring at its northern and southern sides, just 
li ke those seen i n NG C 7552. 

iKohno et al.l ('1999') reported observations of the cir- 
cumnuclear starburst ring of NGC 6951 in both ^^CO (1 

- 0) and HCN (1 - 0) with resolution of -500 pc in both 
lines. They found twin CO knots at or close to the loca- 
tion where the dust lanes meet the ring, and twin HCN 
knots displaced downstream (assuming that the spiral 
arms are tra iling) from the CO knots. On this basis, 
iKohno et al.l ([1999) suggested that the molecular gas at 
the locations where the dust lanes meet the ring — the 
twin CO knots, comprising gas channeled inwards by the 
bar — is too turbulent to form dense molecular clouds. 
Instead, it is only when this gas moves from xi orbits 
(closely traced by the dust lanes) beyond the oILR to 
X2 orbits within the oILR (the circumnuclear starburst 
ring) that dense molecular gas — seen as the twin HCN 
knots — are able to form; the velocity dispersion along 
X2 orbits is predicted to be much smaller than that along 
xi orbits, and hence the molecular gas along X2 orbits is 
much more susceptible to gravitational instabilities (see 
their Fig. 17). This scenario is not supported in our ob- 
servation at least in a global sense as we find that the 
molecular gas in the circumnuclear starburst ring is very 
unsta ble to co l lapse. 

Krip s etldl (120071 ) have observed the circumnuclear 
starburst ring of NGC 6951 at significantly higher angu- 
lar resolution (—240 pc) and sensitivity also in ^^CO (2 

- 1) and HCN (1 - 0). An examination of their results 
(see their Fig. l) reveals no discernible displacement be- 
tween the centroids of the CO knot and dominant HCN 
knot in the northern part of the ring, and at best a small 
displacement between the centroids of the CO knot and 
dominant HCN knot in the southern part of the ring. 
(There is a weaker HCN knot located just downstream 
of each dominant HCN knot: together, these knots, when 
convolved to the lower angular resolution attained in the 
observation by .Kohno et al.. (1999.) . produce an apparent 



Dense Molecular Gas and Stars at Circumnuclear Starburst Ring 



11 



angular displacement between the HCN and CO knots.) 
The centroids of all these knots are located at or close 
to where the dust lanes meet the ring, just as we see 
for the two HCN knots in the circumnuclear starburst 
ring of NGC 7552. It is possible that collisions between 
gas clouds at the intersection of xi orbits (closely traced 
by the dust lanes) and the X2 orbits (the circumnuclear 
starb urst ring) give rise to re latively dense molecular gas 
(e.g.. ICombes fc GerinI IT985I ) . although the relative col- 
lision velocity must not be so high (< 50 km s~^) as to 
dissociate or ionize the molecular gas. Alternatively, or in 
addition, turbulen ce generated by shocks help compress 
the molecular gas (|Elmegreenlll994f ) ; turbulence can help 
support giant molecular clouds against collapse on large 
scales, and promote the formation of dense condensations 

— from which stars form — on small scales. 

As sketched in Figu re [T7| both our o bservations of 
NGC 7552 and that by lKrips et all (|2QQ7[ ) of NGC 6951 
therefore su ggest an important revision to the picture 
proposed bv lKohno et al.l (jTUP). In both these galaxies, 
we suggest that colHsions between diffuse (CO) molec- 
ular gas being channeled inwards along the dust lanes, 
turbulence generated by shocks, along with gravitational 
instabilities perhaps help promote the formation of dense 
molecular gas at or close to the intersection of the xi 
and X2 orbits where the dust lanes meet the circumnu- 
clear starburst ring. Indeed, for NGC 7552, the ratio in 
brightness temperature of the HCN (1 - 0) to ^^CO (2 

- 1) emission averaged over the ring is ~0.18, whereas 
that at the stronger ^^CO (2-1) western arm has a 3a 
upper limit of ~0.12. This indicates that the formation 
of dense molecular gas is promoted where the dust lanes 
meet the ring, relative to regions immediately beyond, 
along the same (at the very least for the western) dust 
lanes. 

Our revised picture also helps explain another observed 
property of the circumnuclear starburst ring in both 
NGC 7552 an NGC 6591, and that is the displacement 
between the centroids of their HCN and radio continuum 
knots. Such a displacement was pointed out above for 
NGC 7552 (see §3.41 and Figure [5j), and is obvious even 
through a curs ory visual co mparison of the molecular line 
maps made bv lKrips et al.i (2007) with the radio contin- 
uum maps made bv'Saiki a et al.l ((2002) in NGC 6951. If 
we assume that the spiral arms in NGC 7552 are trail- 
ing as assumed bv iKohno et al.l ([1999) for NGC 6951, 
then the radio continuum knots are located downstream 
of the HCN knots in both NGC 7552 and NGC 6951. 
The most straightforward explanation for this displace- 
ment is that it reflects the timescale between the for- 
mation of massive stars from dense molecular gas at or 
close to the intersection of the xi and X2 orbits, and 
the eventual demise of these stars in supernova explo- 
sions. This argument can be checked by calculating the 
time it would take for the rotation of the circumnuclear 
starburst ring to have carried the radio continuum knots 
downstream away from the present location where the 
dust lanes meet the ring (specifically, the present loca- 
tion of the HCN knots), keeping in mind that during this 
time the pattern speed of the bar also carries the dust 
lane downstream with respect to the ring, albeit by a 
smaller angular distance. 

Figure [18] shows the 3-cm radio continuum map of Fig- 



ure[3]and HCN integrated intensity map of Figure [5] plot- 
ted in polar coordinates. The northern radio and HCN 
knots have a deprojected angular offset 0'.'47 ± 0'.'05 (~50 
pc), and the southern radio and HCN knots a deproje cted 
angular offset of 0'.'67±0'.'08 (-70 pc). As shown in ^ET| 
the average rotational velocity of both sides of the cir- 
cumnuclear starburst ring is about 218 km s~^, and the 
pattern speed of the bar about 38 km s~^ kpc~^, corre- 
sponding to 19 km s~^ at the ring (0.5 kpc); their net 
velocity is therefore about 199 km s~^. The propaga- 
tion time between the northern radio and HCN knots is 
therefore (2.4± 0.3) x 10^ year, and that between the 
southern radio and HCN knots (3.4 ± 0.4) x 10^ year. 
These timescales are about an order of magnitude shorter 
than the 5-10 Myr timescale for OB stars to explode as 
supernovae. It is possible that the true pattern speed is 
significantly higher and/or the rotational velocity of the 
ring significantly slower than what we inferred, although 
it seems unlikely that both these speeds have been un- 
derestimated by up to an order of magnitude. 

In the case of NGC 6951, the centroids of the twin ra- 
dio continuum knots imaged by Saikia et al. (2002) coin- 
cide closely with the cen t roids of the twin HCN knots as 
imaged by iKohno et al.l ()1999f ) , who estimated a propa- 
gation time between their observed twin CO and HCN 
knots of a few million years. This is in better agreement 
with, although as in NGC 7552 on the shorter side of, the 
timescale for the demise of massive stars with a Salpeter- 
like IMF. We point out that the model proposed above for 
the formation of dense gas and stars in the circumnuclear 
starburst ring can be more cleanly applied to NGC 6951 
than NGC 7552. The latter exhibits not just two but five 
radio continuum knots in its circumnuclear starburst ring 
(see Fig [31 right panel), indicating that the actual pat- 
tern for star formation — and by extension formation 
of dense molecular gas — is more complicated than can 
be fully captured by the simple model proposed in Fig- 
ure [iTl In reality, dense molecular gas may form in the 
dust lanes, as they have been found in some bar g alaxies 
with HCN (1 - 0) observations (e.g., iLeon et al.(l2008l), 
before being channeled into the circumnuclear starburst 
ring (albeit, as shown in the case of NGC 7552, not as 
efficiently as at the intersection of these dust lanes with 
the ring), and continue to form in the circumnuclear star- 
burst ring downstream from where the dust lanes meet 
the ring (e.g., through gravitational instabilities). 

Even in the simpler case of NGC 6591, let alone 
NGC 7552, there are further complications that need to 
be taken into account. The formation of dense molec- 
ular gas at any given location may not proceed at a 
constant rate but is modulated in time. Some star for- 
mation may therefore occur even before gas is channeled 
into the ring and continue further downstream away from 
where the dust lanes meet the ring. In addition, we are 
not able to separate thermal and non-thermal contribu- 
tions to the radio continuum at this stage, and thus can- 
not attribute all the emission to a post-supernova stage. 
Condon (1992) suggests a simple way to separate the 
thermal emission from the total observed radio flux at 
a given frequency based on the spect ral i ndex. For a 
spectral index of -0.89 derived in the §3.3( the fraction 
of thermal emission is about 40% at 3 cm (—10 GHz). 
If thermal (i.e., free- free emission from H II regions) is 
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enhanced at or just downstream of the HCN knots due 
to star formation, then the convolution (due to the lim- 
ited angular resolution of our observation) of this emis- 
sion with non-thermal emission produced further down- 
stream by supernovae would tend to shift the centroid 
of the combined emission upstream towards the HCN 
knots. The result is a smaller apparent propagation time 
between the HCN and radio continuum knots than is ac- 
tually the case. On the other hand, the bright millimeter 
lines seen at the con tact poin ts may reflect high veloc- 
ity dispersions (Down es"et al. l 1993) and therefore large 
line widths, and the concentrations of dense molecular gas 
actually peak somewhat downstream of the HCN knots. 
In this case, the result is a larger apparent propagation 
time between the HCN and radio continuum knots than 
is actually the case. Finally, the IMF for star formation 
in the ring may be top-heavy, resulting in a shorter than 
usual timescale between the formation and death of the 
majority of massive stars (specifically, the time at which 
the radio intensity of all their supernovae combined peaks 
following the formation of these sta rs at the same epoch ) 
compared with a Salpeter-like IMF. Papadopoulos (2010) 
and Papadopoulos et al. (2011) have pointed out that, 
in regions of high star- formation densities, the result- 
ing high cosmic-ray densities can heat molecular gas to 
much higher temperatures than is the case in Galactic 
GMCs, in agreement with the relatively high tempera- 
tures than we infer for the molecular gas in the starburst 
ring of NGC 7552. The higher temperature of the molec- 
ular gas requires a correspondingly higher characteristic 
Jeans mass for star formation in this gas, resulting in: 
(i) stars with larger masses than are found in Galactic 
GMCs; and (ii) a top -heavy IMF as low-mass star for- 
mation is suppressed (Papadopoulos e t "ani2011f ). 

Our revised model of the sequence of star forma- 
tion in this circumnuclear ring is consisten t with the 
pearls on a string scenario of BokerelaD (j2008f ) (see 
their Figure 7), which is based on their near infrared 
(H- and K-band) observations along with the age identi- 
fication using the intensity ratios of emission lines. The 
pearls on a string model proposes that the overdense 
region is close to where the gas enters the ring, and 
therefore where the starburst is triggered. The star clus- 
ters formed subsequently will continue their orbit along 
the circumnuclear ring. In this scenario, a bi-polar age 
gradient is expected, as has been seen in some galax- 
ies, based on infra red an d optical st udies, such as NGC 
1068 (Da vies et al.1ll998h NGC 7771 (ISmith et al.i ri999). 
MlOO (.Rvder et al.l 120011 lAllard et a ljl2006f). IC 4933 
(jRvder et al.ll2010f ). and M83 (KnaDen 'eraIll2010f ). 

The infrared age dating method has been applied to 
NGC 7552 by Brandl et al. (2012), who identify nine un- 
resolved mid-infrared peaks (star clusters) in the star- 
burst ring. There are tentative signs of an age gra- 
dient, in that the relatively young clusters are located 
near the north contact point, while the older clusters are 
found more than a quarter of the way around the ring 
away from the contact points. This result is qualita- 
tively consistent with the age gradient expected by the 
pearls on a string scenario and our revised model, al- 
though it must be noted that the age difference among 
the clusters is rather small at between o nly 5.5 and 6.3 
Myr (see Figure 10 in lBrandl eraIll2Q12f ). 

Even though the aforementioned bi-polar age gradi- 



ents are fo u nd in the circumnuclear r ing s urveys of 
iBoker et al.l (j2008D and iMazzuca etld] (j2008f ). the de- 
tection rate of such age gradient patterns is only 60% 
and 43%, respectively, in these two surveys. We can- 
not therefore rule out the possibility of other scenarios 
for star formation, in addition to that described above, 
in the ci rcumnuclear s tarbu rst ring of NGC 7552. For 
example, iBoker et al.l ()2008l ) suggest a second Popcorn 
model in which the gas clo uds become g;rav itationally 
unstable due to turbulence (Elmegreen" "1994) , and col- 
lapse at more random times and locations. Accordingly, 
there is no clear a ^e gradient expect ed in the Popcorn 
model, as found bv lHsieh et al.l ()2011f ) for the circumnu- 
clear starburst ring of NGC 1097. 

In NGC 7552, we detect five radio continuum knots, 
not just two. In NGC 6591, apart from the two dom- 
inant HCN(1 - 0) peaks, there are two other HCN(1 - 
0) knots in the circumnuclear starburst ring. Thus, even 
for NGC 7552 and NGC 6591, the overah picture for the 
formation of dense molecular gas and stars in circumnu- 
clear starburst rings cannot be as simple as that sketched 
in Figure [ITl In reality, the formation of dense molecu- 
lar gas may be promoted at the contact points, but also 
elsewhere in the circumnuclear starburst ring and star 
formation also occur there. Both the pearls on a string 
and Popcorn models may operate simultaneously in any 
given circumnuclear starburst ring, complicating any sin- 
gle interpretation. What we present here is a zeroth- 
order model; the true picture requires observations at 
higher angular resolutions to properly study all locations 
where dense molecular gas is forming. 

7. SUMMARY AND CONCLUSIONS 

NGC 7552 is a barred gal axy th at, like a significant 
fraction (~20%-50%; Kna penll2005l ) of such galaxies, ex- 
hibits a circumnuclear starburst ring. A pair of dust 
lanes can be seen running along the bar on opposite 
sides of the center of the galaxy, and at their inner ends 
connect with the northern and southern parts of the 
circumnuclear starburst ring. The circumnuclear ring 
currently contributes 50% - 80% of the galaxy's entire 
star formation, and has a gas consumption timescale 
intermediate between normal and extreme star-forming 
(ultraluminous-infrared) galaxies. 

To better understand the sequence and circumstances 
under which dense molecular gas and stars form in this 
ring, we have observed the center of NGC 7552 in molec- 
ular lines tracing relatively diffuse molecular gas (^^CO 
and ^^CO J = 2 - 1), which is the characteristic of Galac- 
tic giant molecular clouds; and relatively dense molecular 
gas (HCN and HCO+ J = 1 - 0), which is the charac- 
teristic of dense condensations in giant molecular clouds 
where stars form. We find that emission from: 

• the dense molecular gas tracers coincide with 
the circumnuclear starburst ring, and exhibit two 
peaks located at or close to where the dust lanes 
meet this ring, and 

• the diffuse molecular gas tracers also coincide with 
the circumnuclear starburst ring, and ^^CO (2-1) 
also with the innermost ends of the two dust lanes. 

Using the LVG approximation, we demonstrated that: 
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• the molecular gas as observed averaged over the cir- 
cumnuclear starburst ring is both relatively dense 
(10^^^ cm~^) and warm (well over 100 K), sug- 
gesting that massive stars have indeed formed from 
this gas, and 

• the diffuse molecular gas in the circumnuclear star- 
burst ring is unstable to gravitational collapse, pro- 
viding at least one pathway for the formation of 
dense molecular gas. 

We recreated maps of the H I emission from NGC 7552 
using the data taken by |Dahlem (2005), but at a higher 
angular resolution than those previously published. We 
used this map, along with our maps in molecular lines, 
to derive an approximate rotation curve for NGC 7552, 
and computed the predicted locations of dynamical res- 
onances under the usual assumption that the outer ends 
of the bar terminate at the CR. We found that, within 
the statistical uncertainties, our results: 

• are consistent with the idea that the OLR coincides 
closely with the outermost reaches of the twin spi- 
ral arms (which emerge from the ends of the bar) 
in NGC 7552, as theory predicts, and 

• show that the circumnuclear starburst ring likely 
lies between the ilLR and oILR, as theory also pre- 
dicts. 

We then examined what our results imply for the for- 
mation of dense molecular gas and stars in the circum- 
nuclear starburst ring of NGC 755 2, specifically in the 
context of the picture proposed by iKohno et al I fl999) 
based on their observations of NGC 6951. The results 
and conclusions are as follows: 

• In both NGC 7552 and NGC 6951, emission from 
the dense molecular-gas tracers peaks at or close to 
the location where the pair of dust lanes that run 
along the bar meet the circumnuclear starburst ring 
(in the case of NGC 695 1 , as l ater improved ob- 
servations by iKrips et al.l ()2007[ ) also show) . This 
result suggests that the formation of dense molecu- 
lar gas is promoted where the xi orbits (as closely 
traced by the dust lanes) and X2 orbits (circum- 
nuclear starburst ring) intersect and where colli- 
sions between gas clouds is therefore more likely, 
rather than further downstreani along the ring as 
proposed by 'Kohno et ajj (| 19991 ) where the veloc- 
ity dispersion of X2 orbits is predicted to be much 
smaller than that of xi orbits and hence the gas 
more susceptible to gravitational instabilities. In- 
deed, we observe that the formation of dense molec- 
ular gas is promoted where the dust lanes meet 
the ring compared with regions immediately be- 
yond along the same dust lanes. 

• Our revised model as sketched in Figure [TTl also 
naturally explains the observed displacement be- 
tween the centroids of the dense molecular gas 
and radio continuum knots in both NGC 7552 and 
NGC 6951. This displacement ought to reflect the 
timescale between the formation of massive stars 
in the dense gas located where the dust lanes meet 



the circumnuclear starburst ring, and the demise of 
these stars as supernovae. In the case of NGC 6951, 
the estimated timescale between the centroids of 
these knots is a few million years, at the shorter 
end of the 5-10 Myr timescale for the majority of 
massive stars formed with a Salpeter IMF to be- 
come supernovae, and in the case of NGC 7552 is 
an order of magnitude shorter still. This may re- 
flect uncertainties in the inferred angular velocities 
of the ring and bar, together with additional com- 
plexities mentioned below or a top-heavy IMF. 

Finally, we noted that the observed properties of the 
circumnuclear starburst ring in NGC 7552 are more com- 
plicated than can be fully explained by the simple model 
proposed earlier and sketched in Figure [TTl This model 
can be more cleanly applied to NGC 6951, whereas the 
circumnuclear starburst ring of NGC 7552 exhibits not 
just two, but five, knots in radio continuum. Our obser- 
vations do not have sufficiently high angular resolution 
to spatially separate these knots in the dense (or diffuse) 
molecular gas tracers, so that the picture that emerges 
from our analyses represents a zeroth-order model for 
the formation of dense gas and stars in the circumnu- 
clear starburst ring of NGC 7552. Observations at higher 
angular resolution and sensitivity are necessary to tease 
out details not incorporated into our model, which itself 
is based on an important revisio n of the model originally 
proposed bv lKohno etHI ^M9t ) for NGC 6951. 
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TABLE 1 

Observational parameters of radio observations. 





H I 


6 cm 


3 cm 


HCN(l-O) 


HCO+(1-0) 


i3CO(2-l) 


i2CO(2-l) 


Rest Frequency (GHz) 


1.420 


4.786 


8.640 


88.631 


89.188 


220.398 


230.538 


Telescope 


ATCA 


ATCA 


ATCA 


ATCA 


ATCA 


SMA 


SMA 


Projected baseline (m) 


245 - 6000 


153 - 6000 


153 - 6000 


30 - 200 


30 - 200 


16 - 69 


16 - 69 


Field of view 


33' 


285'' 


285" 


36" 


36" 


55" 


55" 


Velocity resolution (km s"-*^) 


14 






15 


15 


18 


10 


Synthesized beam (arcsec^) 


20x20 


2.3x1.3 


1.1x1.1 


2.6x2.0 


2.6x2.0 


6.9x2.8 


7.0x2.8 


P.A. (deg) 











80 


80 


-8.9 


-11.9 


Reprocessed data (yes/no) 


y 


y 


y 


n 


n 


n 


n 


Reference 


1 


2 


2 











^Dahleml (|200^) 
^Forbes et all (|1994al) 



Database (NED) which is operated by the Jet Propul- 
sion Laboratory, Cahfornia Institute of Technology, un- 
der contract with the National Aeronautics and Space 



Administration. 

Facilities: ATCA, SMA, Gemini. 
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Fig. 1. — Optical r'- (upper left), g'- (upper right), i^-band (lower left) and Ho; (lower right) images of NGC 7552 taken with CMOS 
on the Gemini South telescope. The Spitzer 5.8 /xm contour at 5.5 MJy sr"-*^ (Kennicutt et al. 2003) is overlaid on the ^f^-band image, 
highlighting the orientation of the dust lanes. Likewise, the 3-cm continuum contours (" ^3.31 see text) showing the circumnuclear starburst 
ring are plotted in the Ho; image, with contour steps of 3, 15, and 35 cr, where 1 cr = 66 fiJy beam"-*^. 
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Fig. 2. — Integrated H I intensity (left panel) and intensity- weighted mean H I velocity (right panel), both in color. In the left panel, 
contours are of the i' image shown in Figure [T] In the right panel, contours correspond to the intensity-weighted mean H I velocity, just 
like the color map, and are plotted in steps of -90, -45, 0, 45, and 90 km s"-*^ with respec t to the syste mic velocity. In each panel, a cross 
is plotted at the center of the galaxy, corresponding to the center of the starburst ring ((Forbes et al .ill 994a ). and the synthesized beam 
shown at the lower left corner. 
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Fig. 3. — Radio continuum images of the circumnuclear starburst ring in NGC 7552 at 6 cm (left panel; both contours and greyscale) 
and 3 cm (right panel; contours) recreated from the data taken by and published in Forbes et al. (1994a). In the right panel, the 3-cm 
(contours) is overlaid on the Ha (color scale) image shown in Figure ^ Straight lines in both panels indicate peaks in radio emission, which 
correspond to star-forming knots. Contour levels of the 6-cm image are at 3, 10, 25, 40, 55, 70, and 85cr, where 1 cr = 72 /iJy beam"-*^, and 
contour levels of the 3-cm image are a t 3, 5, 15, 25, 35, and 40 cr, where 1 cr = 66 /iJy beam"-*^. A cross in each panel indicates the center 
of the ring as defined by Forbes et J] (ll994al V The synthesized beam is 2^'0 x 1^'3 at 6 cm and l^^l x 1^1 at 3 cm, and shown at the lower 
left corner of the respective panels. 
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Fig. 4. — Channel maps of the central region of NGC 7552 in HCN (1 - 0) (black contours) and HCO"^ (1 - 0) (cyan contours). Contour 
levels are plotted at -2, 2, 3, 5, 7, and 10 x 4.2 mJy beam"-*^ (rms noise level) for HCN (1 - 0) and x4.5 mJy beam"-*^ (rms noise level) for 
HCO"^ (1 - 0). The velocity of each channel, in units of km s"-*^, is shown in the upper left corner of each panel. A cross is plotted at the 
center of the galaxy in each panel. The synthesized beam is 2(^6 x 2^0 at a position angle of 80°, and is shown at the lower left corner of 
each panel. 
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Fig. 5. — Integrated HCN (1 - 0) intensity (left panel; contours) and intensity-weighted mean HCN (1 - 0) velocity (right panel; contours 
and color scale) derived from the channel maps shown in Figure [U In the left panel, contours indicating the integrated HCN (1-0) 
intensity are plotted at 3, 6, 10, 12, and 15 cr, where 1 a = 0.3 Jy beam~^ km s~^. Greyscale corresponds to the 3-cm continuum image 
shown in Figure [S] but convolved to the same angular resolution as the HCN (1 - 0) image. In the right panel, contours correspond to the 
intensity- weighted mean HCN (1 - 0) velocity, just like the color map, and are plotted in steps of 20 km s"-*^ starting from —60 km s"-*^ 
and ending at 60 km s"-*^ with respect to the systemic velocity. A cross is plotted at the center of the galaxy, and the synthesized beam 
shown at the lower left corner, of each panel. 
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Fig. 6. — ntegrated HCO"^ (1 - 0) intensity (left panel; contours) and intensity-weighted mean HCO+ (1 - 0) velocity (right panel; 
contours and color scale) derived from the channel maps shown in Figure [H In the left panel, contours indicating the integrated HCO+ (1 
- 0) intensity are plotted at 3, 6, 10, 12, and 15 cr, where la = 0.3 Jy beam~^ km s"-*^. Greyscale corresponds to the 3-cm continuum image 
shown in Figure [3] but convolved to the same angular resolution as the HCO+ (1-0) image. In the right panel, contours correspond to the 
intensity- weighted mean HCO+ (1 - 0) velocity, just like the color map, and are plotted in steps of 20 km s"-*^ starting from —60 km s"-*^ 
and ending at 60 km with respect to the systemic velocity. A cross is plotted at the center of the galaxy, and the synthesized beam 
shown at the lower left corner, of each panel. 
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Fig. 7. — Channel maps of the central region of NGC 7552 in ^^CO (2 - 1). Contour levels are plotted at -2, 2, 3, 5, 7, 10, and 
12 X 39 mjy beam"-*^ (rms noise level). The velocity of each channel, in units of km s~^, is shown in the upper left corner of each panel. 
A cross is plotted at the center of the galaxy in each panel. The synthesized beam is x 2^8 at a position angle of -8?1, and is shown 
at the lower left corner of each panel. 
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Fig. 8. — Integrated -"^^CO (2 - 1) intensity (left panel; contours) and intensity- weighted mean ■'^'^CO (2 - 1) velocity (right panel; contours 
and color scale) derived from the channel maps shown in Figure[7l In the left panel, contours indicating the integrated -"^^CO (2 - 1) intensity 
are plotted at 3, 10, 14, and 18cr, where la = 2.3 Jy beam"-*^ km s"-*^. Grayscale corresponds to the integrated HCN (1 - 0) intensity image 
shown in Figure O but convolved to the same angular resolution as the -"^^CO (2 - 1) image. In the right panel, contours correspond to the 
intensity- weighted mean -"^^CO (2 - 1) velocity, just like the color map, and are plotted in steps of 20 km s"-*^ starting from —60 km s"-*^ 
and ending at 60 km s"-*^ with respect to the systemic velocity. A cross is plotted at the center of the galaxy, and the synthesized beam 
shown at the lower left corner, of each panel 
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Fig. 9.— Position-velocity (P-V) diagrams of HCN (1 - 0) (upper left panel), HCO+ (1 - 0) (upper right panel), ^^CO (2 - 1) (lower 
left panel) and ^^CO (2 - 1) (lower right panel) along a position angle of 110°, which corresponds to the orientation of the kinematic 
major axis. Negative positions are located eastwards and positive positions located westwards of the galaxy center. Dashed horizontal lines 
indicate km s~^ relative to the systemic velocity of 1580 km s~^. Contour levels of HCN (1-0) and HCO"'" (1 - 0) are at —2, 2, 3, 5, 7, 
and 10 where 1 a is 4.2 and 4.5 mjy beam-\ respectively, for HCN (1 - 0) and HCO+ (1 - 0). Contours levels of ^^CO (2 - 1) are at 
-2, 2, 3, 5, 7, 10, 12, and 15 a, where 1 a = 39 mJy beam"!, and ^^CO (2 - 1) at -2, 2, 3, 5, 7, 15, 30, and 40 a, where 1 a = 78 mJy 
beam~^. 
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Fig. 10.— Channel maps of the central region of NGC 7552 in ^^CO (2 - 1). Contour levels are plotted at -2, 2, 5, 10, 25, 40, and 
55 X 78 mjy beam"-*^ (rms noise level). The velocity of each channel, in units of km s~^, is shown in the upper left corner of each panel. 
A cross is plotted at the center of the galaxy in each panel. The synthesized beam is 7'.'0 x 2(^8 at a position angle of -11? 9, and is shown 
at the lower left corner of each panel. 
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Fig. 11. — Upper left: Integrated -"^^CO (2 - 1) intensity map in white/black contours plotted at 2, 5, 10, 15, 20, 25, 45, 65, and 
85 X 4.5 Jy beam~^ km s~^ overlaid on the optical i'-band image in grayscale shown in Figure [T] Upper right: Contours of the integrated 
HCN (1 - 0) intensity map (contour levels are the same as those shown in Figure [5]) overlaid on the greyscale i^-band image. Lower left: 
Integrated -"^^CO (2 - 1) (green contours) and HCN (1 - 0) (black contours) intensity maps both overlaid on the i^-band image (grayscale). 
Lower right: Intensity- weighted mean -"^^CO (2 - 1) velocity in color scale, as well as white contours plotted in steps of 20 km s"-*^ starting 
from —60 km and ending at 60 km with respect to the systemic velocity of 1586 km s~^. In each panel, cyan contour is the outline 
of the 5.8 yum emission at 5.5 MJy sr"-*^ from the Spitzer SINGS project published in Kennicut t et al.i ((2003|); this emission traces the dust 
lane visible in silhouette in the optical image. A cross is plotted at the center of the galaxy, and the synthesized beam shown at the lower 
left corner, of each panel. 



26 



Pan et al. 




10"' 10"' 10"' 10"' 10"' 10"' 10"' 10"' 10"" 
Z(''CO)/(dv/dr) Z('^CO)/(dv/dr) 

Fig. 12. — Derived physical properties of the molecular hydrogen gas in the circumnuclear starburst ring based on our LVG analysis. In 
all panels, density of molecular hydrogen gas is plotted as a function of the product Z {^^ CO) / {dv / dr) ^ which for the present circumstances 
has an adopted value of ~4 x 10~^ (vertical line in each panel; see text). Black contours in each panel indicate the -"^^CO (2 - 1) opacity. 
Red contours are the ratio in brightness temperature of HCN (1 - 0) to -"^^CO (2 - 1), which averaged over the ring has a value of 2.0 ± 0.2 
(thick red contour in each panel where visible). The upper left panel are solutions for a kinetic temperature of 10 K; upper right panel 50 K; 
lower left panel 100 K; and lower right panel 1000 K. As can be seen, for the measured line ratio and adopted product Z {^^ CO) / {dv / dr) ^ 
solutions can only be found for kinetic temperatures well above 100 K. 
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Fig. 13. HCN (1 - 0)/^'^CO (1 - 0) versus HCN (1 - 0)/HCO+ (1 - 0) integrated intensity ratios of pure Seyfert and starburst galaxies. 
Solid circles indicate Seyfert galaxies, open squares starburst gal axies, and an open star the position of NGC 7552 (starburst) in this plot. 
Data other than for NGC 7552 are taken from .Kohng (,2005 ), Kohno et al.l (,2007 , 2008) and ?. 
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Fig. 14. — Rotation curve, projected for an inclination of 25°, shown in red curves overlaid on a P-V diagram derived from our ^^CO (2 
- 1) (—1 to +1 kpc) and H I maps —1 and ^ 1 kpc) at a position angle of 110° shown in contours. Shaded regions indicate locations 
where the two tracers connect, and at which the rotation curve fail to track (see text). The typical projected uncertainty of the rotation 
curve obtained by the envelope tracing method is ±15 km s"-*^, and illustrated at the end of the approaching side. Contours in -"^-^CO (2 - 
1) are plotted in steps of -2, 2, 3, 5, 10, 20, 30, 40, and 45 cr, where 1 cr = 78 mJy beam"-*^, and contours in H I diagram in steps of —2, 2, 
3, 5, 7, 10, 15, and 20cr, where 1 cr = 1.1 mJy beam~^. 
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Fig. 15. — Deprojected rotation curves derived for NGC 7552 from our ^^CO (2 - 1) maps (< 1 kpc) and H I maps (> 1 kpc). Left 
panel shows the rotation curve within 1 kpc, and right panel the entire rotation curve up to 10 kpc. Solid and open black circles denote 
the rotation curve of the receding and approaching sides, respectively. We assumed a constant inclination of 25° and position angle of 110° 
for the kinematic major axis . Red circles show the average rotation curve at the radii where both the receding and approaching sides can 
be measured. The uncertainties of the individual and average rotation curves are indicated in the lower right corner of each panel. 



30 



Pan et al. 



Radius [arcsec] 

18.7 37.4 56.1 74.7 




2 4 6 8 

Radius [kpc] 

Fig. 16. — Diagram used to infer the locations of dynamical reso nan ces. Black solid circles indicate the angular velocity of the galaxy, 
based on an interpolation of the rotation curve shown in Figure fTSl Individual measurements that define the fl — k/2 (ILR) curve are 
plotted as triangles, the fl — k/A (UHR) curve as diamonds, and the k/2 (OLR) curve as squares, where k is the epicyclic frequency of 
the orbit. Resonances occur where the values of these curves are equal to the pattern speed of the bar, indicated by the thick horizontal 
line. 
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Fig. 17. — A sketch showing our zeroth-order picture of the formation of dense gas (HCN peaks) and synchrotron emission (3-cm peaks) 
resulting from supernova explosions of stars that formed from the dense gas in the circumnuclear starburst ring of NGC 7552. This sketch 
is based on a revision of the picture proposed by Kohno et al. (1999) for NGC 6591. 
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Fig. 18. — Integrated HCN (1 - 0) intensity of Figure [5] (upper panel) and 3-cm continuum intensity of Figure [3] (lower panel) plotted in 
polar coordinates. The abscissa is polar angle measured from east, and the ordinate radius from center. Crosses in both panels indicate 
the centroids of prominent pairs of knots in the emission. 



